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ABSTRACT 

In this letter we report on turbulent acceleration of the dissipation of magnetic field in the postshock re- 
gion of a Poynting flux-dominated flow, such as the Crab pulsar wind nebula. We have performed two- 
dimensional resistive relativistic magnetohydrodynamics simulations of subsonic turbulence driven by the 
Richtmyer-Meshkov instability at the shock fronts of the Poynting flux-dominated flows in pulsar winds. We 
find that turbulence stretches current sheets which substantially enhances the dissipation of magnetic field, and 
that most of the initial magnetic field energy is dissipated within a few eddy-turnover times. We also develop a 
simple analytical model for turbulent dissipation of magnetic field that agrees well with our simulations. The 
analytical model indicates that the dissipation rate does not depend on resistivity even in the small resistivity 
limit. Our findings can possibly alleviate the cr-problem in the Crab pulsar wind nebulae. 
Subject headings: magnetic fields, magnetohydrodynamics (MHD), relativistic processes, shock waves, plas- 
mas, turbulence 



1. INTRODUCTION 

Magnetic fields play an important role in various as- 
trophysical phenomena. In particular, Poynting flux- 
dominated plasmas have been studied extensively as ori- 
gins of high energ y astrophysical phenom ena, e.g., ul- 
tra relativistic j ets (Meszaros & Rees 1 1 19971) . gamma ray 
bursts (GRB) dLvutikov & Blandford I 120031) . and pulsar 
winds dKennel & Coroniti 1 1 1 984allbb . The key point of these 
models is that most of the energy is transferred as Poynting 
flux from the central engine first, and then converted into 
the thermal energy later. This makes it easy to form high 
Lorentz factor conditions and well collimated jets since the 
inertia of the plasma is negligible. However, it is widely 
known that it is very difficult to sufficiently dissipate the elec- 
tromagnetic energy by simple collisional Ohmic dissipation 
within the observationally indicated characteristic times of 
phenome na, so that many alternative mechanism s have been 
proposed (Meie rl2003HGiannios & Spruit 120 06). Analogous 
problems can also be found in non-rel ativistic phenomena, 
such as solar flares (Sturrock et al. 1984), and it seems to be 
a generic problem in high magnetic Reynolds number media. 
Therefore, research into mechanisms of efficient dissipation 
of electromagnetic energy is extensive. 

The cr-problem is a puzzle of energetics be- 
tween t he Crab pulsar wind a nd pulsar wind nebula 
(PWN) dKennel & Coroniti I[l984allbh . Theoretical studies of 
the Crab pulsar magnetosphere suggest that the pulsar wind is 
a Poynting flux-dominated flow with a ~ 10 4 , where a is the 
ratio of the Poynting energy flux to the particle energy flux. 
Furthermore, if one assumes an axisymmetric stationary rel- 
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ativistic magnetohydrodynamic (RRMHD) flow, the plasma 
continues to be Poynting flux-dominated far beyond the fast 
magnetosonic point. On the other hand, observations of the 
expansion speed and luminosity of the Crab PWN indicate 
that the plasma should be particle energy -dominated near the 
inner e dge of the nebula as a ~ 3 x 10~ 3 dKennel & Coronitil 
Il984alfb1) . To resolve this problem, efficient dissipation 
processes are necessary in the wind region. One of the we ll- 
known models is the "striped wind model" dMichel II 19711) in 
which oblique rotating pulsars have many radially outgoing 
current sheets around the equatorial plane in the wind region. 
This model was expected to resolve the cr-problem through 
dissipation in these current sheets. It was, however, shown 
that the magnetic field energy cannot be fully dissipated and 
Poynti ng energy still dominates at the inner edge of the Crab 
PWN dLvubarskv & Kirkll200lt lkirk & Skiaeraasenll2003h . 
Concerning this problem, Petri & Lyubarsky (2007) and 
Sironi & Spitkovsky (2011) have performed particle-in-cell 
(PIC) simulations that shown the magnetic field dissipation 
rate is substantially enhanced by the interaction of a weak 
precursor MHD shock wave and a current sheet that induces a 
driven magnetic reconnection. These works largely advanced 
our understanding of the cr-problem, which focused mainly 
on kinetic effects. In this letter, we show that hydrodynamical 
effects can also provide an additional mechanism of magnetic 
field dissipation behind the precursor MHD shock wave. The 
preshock region of the pulsar wind can be inhomogeneous 
due to many possible sources of fluctuation, such as magnetic 
reconnection at the Y-point close to the light cylinder; the 
fluctuations in the preshock flow can trigger instabilities 
and generate turbulence in the postshock region, such as the 
Richtmyer-Meshkov instability (RMI) driven by the shock- 
density fluctuation interaction (Giacalone & Jokipii 2007; 
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Sironi & Goodman 2007; Inoue et al. 2009,2010,201 1,2012). 
In such a turbulent medium, current sheets in the striped wind 
can be stretched by turbulent vortical motions, which can 
enhance the dissipation of magnetic energy. In the following, 
we consider decaying turbulence by driving it only in the 
initial setup, as in the case of shock-induced turbulence by 
the RMI in which it is only induced at the shock front. Such 
mechanism of magnetic field dissipation is analogous to the 
turbulent reconnection proposed by Lazarian & Vishniac 
1999 (see also Kowal et al. 2009). While turbulent reconnec- 
tion assumes externally driven turbulence, our results show 
that the fast dissipation of magnetic fields proceeds even 
without external forcing, which strengthens the applicability 
of our mechanism. 

2. ENERGETICS TO DRIVE TURBULENCE 

In the aforementioned scenario, the turbulence induced in 
the postshock region is driven by thermal and magnetic en- 
ergy released in the high-cr shock. Hence, the velocity dis- 
persion of the turbulence is expected to be less than the slow 
magnetosonic velocity behind the high-<7 shock wave. In this 
section, we estimate the upper limit of the velocity dispersion 
of the turbulence. 

In the case of the high-cr perpendicular RMHD shock, 
the physical variables describing the pos tshock region in the 
shock rest frame are written as follows dKennel & Coronitil 
fl984ah : 

B 2 ^B U pzlz — PHi, 72^%/cr, k B T 2 /mc 2 ^^i/Si/a, 

(1) 

where subscripts 1 and 2 refer to preshock and postshock, 
respectively, p is the mass density in the plasma comoving 
frame, 7 is the Lorentz factor in the shock front rest frame, 
m and c are the particle rest mass and light velocity, respec- 
tively, and cr is the magnetization parameter in the preshock 
flow defined as: 

o = B\/^p l c 1 1 l (2) 

From the above relations, the thermal energy density in the 
postshock rest frame is given by 

E th = Pl c 2 tf/Sa(r-l), (3) 

where T is the adiabatic index, and here we use the value in 
the relativistic limit, T = 4/3. In the postshock region of the 
high-cr flow, the dominant inertia of the plasma is due to the 
magnetic field energy. Thus, the kinetic energy density of the 
turbulence in the postshock rest frame can be evaluated as 

Eturb = [P2C 2 + Bo/ 4Tr"fl]^, urb ("f,urb- 1) — (#2 Av, H rf,) 2 /87r<T, 

(4) 

where Av tur b is the velocity dispersion of the turbulence. 
Since the kinetic energy density of the turbulence should be 
less than the thermal energy density {E tur j, < £,/,), the maxi- 
mum velocity dispersion of turbulence induced by the shock 
can be evaluated as 

Av turb /c s < (5) 

where we have used a = B 1 l /ATtp\c 2 y L , and the relativistic 
sound speed c s = c/^/3. 

According to t he oretic al studies of striped wind s 
dLvubarskv & Kirk I l200lt iKirk & Ski aeraasen ] E003), 
the magnetization parameter immediately upstream of the ter- 
mination shock is estimated as a ~ 100. In this case, Eq. (0 
suggests that turbulence with Av rur b < 0.1 c s can be induced 



in the postshock flow. Note that, according to the relativistic 
MHD simulations by Inoue et al. (201 1), the kinetic energy of 
turbulence can be comparable to the released particle kinetic 
energy when the preshock density dispersion is comparable 
to the mean density, and the velocity dispersion of turbulence 
depends linearly on the density dispersion. This suggests that 
the maximum velocity dispersion given in Eq. (5) is the case 
when the preshock density dispersion is comparable to the 
mean density. 

In the above discussion, we consider the RMI as a candidate 
mechanism for driving turbulence in the postshock region. 
There are, however, other instabilities in the pulsar wind, for 
example, a Rayleigh-Taylor or Kruskal-Sch warzschild type 
instability of current sheets (Lvubarskv l l2010l) that may also 
drive turbulence in a striped wind. 

3. NUMERICAL SETUP 

We model the evolution of the current sheet in turbulence 
using 2-dimensional resistive relativistic magnetohydrody- 
namics (RRMHD). Since this calculation includes relativis- 
tically hot plasma, turbulence and magnetic dissipation, we 
need a highly accurate numerical scheme that can treat rela- 
tivistic Ohmic dissipation and the evolution of magnetic field 
in a turbulent medium. We use a multi-dimensional exten- 
sion of the resistive RMHD scheme developed by Takamoto 
& Inoue (201 1). 1 . This scheme uses appropriate character- 
istic velocities for calculating numerical fluxes, i.e., fast mag- 
netosonic velocity and Alfven velocity, so that it describes 
turbulent flows accurately (see Takamoto & Inoue 2011 for 
detail). We calculate the RRMHD equations in a conser- 
vative fashion, and the mass density, momentum, and en- 
ergy are also conserved within machine round-off error. For 
the equation of state, we assume a relativistic ideal gas with 
h = 1 + (r/(r- l))(p gas /p) where T = 4/3 and h is the relativis- 
tic enthalpy. The resistivity, 77, is assumed to be constant. We 
examine the following five cases: r//c 2 t = 1,2,4,8, 16 x 10" 4 
to study the effects of resistivity, where fo is the light crossing 
time for the numerical domain. 

In our simulations, we model turbulent flow in the post- 
shock region generated by a Poynting flux-dominated flow 
in the post-shock rest frame. Since we consider the situ- 
ation of a very strong magnetic field, turbulence is highly 
sub-Alfvenic. Recent studies have shown that sub-Alfvenic 
turbulence is highly anisotropic and that the structure of 
turbul ent eddies is elongated along the mean magnetic 
field dGoldreich & Sridharll99l iMaron & Goldreich lEOOll 
Cho et al. I 120021: llnoue et al. I 1201 lb IZrake & MacFadyen I 
20121) . In the limit of a strong magnetic field, turbulence 
can be approximated as two-dimensional (2D) in the plane 
perpendicular to the mean magnetic field. Hence, as a first 
step in this kind of simulation, we perform two-dimensional 
simulations in the plane perpendicular to the magnetic field. 
This restriction greatly reduces the computational cost and en- 
ables us to study the dependence of basic parameters, such 
as the resistivity and the strength of turbulence, on mag- 
netic field dissipation. The two dimensional approximation 
also limits magnetic energy dissipation to Ohmic dissipation, 

1 It is widely known that we must take into account the evolution of the 
electric field for satisfying causality; this essentially changes the mathemati- 
cal properties of the governing equation s, and makes it difficult to develop an 
accurate numerical scheme of RR MHD (Komissarov 2007; Palenzue la et al. I 
2009; Takamoto & Inoue 2011). An analogous problem also exists in the 
description of other dissipat i on mechanisms in the framework of relativity 
Miscock & Lindblomll9 85; Takamoto & Inutsuka 2010, 20TQ). 
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so that neit her Sweet-Parker nor Petschek typ e magnetic re- 
connection (Biskamp 20001 iPriest & Forbes I [2000) can oc- 
cur. However, the characteristic timescales of the relativis- 
tic tearing instability and the Sweet-Parker reconnection is 
given by r rec ~ (L/c J 4)v / Rm, where c A is the Alfven veloc- 
ity, Rm is the magnetic Reynolds number defined as Rm 
= Lv/rj, L is the characteristic length sc ale, v is the char - 
acteristic velocity and 77 is the resistivity ( Lyuba rsky] 120051) . 
On the other hand, the timescale of turbulence is given by 
eddy-turnover time L/Av, where Av is the velocity dis- 
persion. The ratio of these two timescales is T lur b/T rec ~ 
(L/Av)/(L/c A \fRm) — (c A / Av)/ \/Rm, indicating that the 
effect of turbulence becomes faster th an t he Sweet-Parker 
type reconnection when Av/c A > 1 / \/Rm 2 . Note that in 
the realistic three-dimensional case, even if the Sweet-Parker 
type reconnection is slow, the turbulence might induce "turbu- 
lent r econnection" (lLazarian & Vishniac |[l99 9; Kowa fet al. I 
2009). The 3D effect will be studied in our future paper. 

For our numerical calculations, we prepare a square do- 
main of [L x x L y ] = [1 x 1], and divide into the homoge- 
neous numerical meshes of N x x N y = 1024 2 . We set periodic 
boundary conditions for the x and y-directions. For the ini- 
tial condition, we consider the following magnetic field pro- 
file: B z = B for x < 0.25L V , or x > 0.15L X , and B z = -B for 

0. 25L A - < x < 0.75L A . In addition, we consider the relativis- 
tically hot postshock region of the Poynting flux-dominated 
cold preshock flow. According to the shock jump conditions 
of high-er flow (see Eq. JT]), the plasma (3 of the postshock 
region is approximately given by (3 ~ er~ 2 where f3 is defined 
as: (3 = p gas /(B 2 /871"). Thus, for a flow of a ~ 100, we have 
(3 ~ 10" 4 . Unfortunately, our code having strict conservation 
properties cannot stably solve such a low (3 plasma. In or- 
der for the stability, we choose (3 = 4, i.e., we adopt a plasma 
with higher thermal pressure. The bottom panel of Fig. |2] 
is the results of simulations in the case of (3 = 0.5, 1,4 along 
with Av = 0.08c s and k c = 3. This shows that small number of 
plasma (3 < 0.5 produces essentially the same dissipation rate, 
because the rate of the dissipation does not essentially depend 
on the f3 (see, eq. [8] below) 3 . Note that the postshock state 
of the high-cr flow is relativistically hot (c s = c/v3), so that 
our choice of a larger plasma (3 does not lead to an overesti- 
mation of the sound speed in the postshock region. 

In this study, we consider so-called "decaying turbulence", 

1. e., the turbulence is driven only initially. This is because the 
shock induced turbulence (e.g. by the RMI) is expected to be 
triggered only at the shock front. We set a divergence-free 
initial turbulent velocity field whose one-dimensional power 
spectrum is flat with a cut off wave number k c = 3/L x or 5/L x , 
i.e., the initial turbulence is driven only at large scales. We 
study the cases of the following five initial velocity disper- 
sions: Av/c s = ^/ly^/c, = 0.05,0.08,0.1,0.3,0.5 where c s 
is the sound velocity. These parameters are motivated by the 
energetics of the turbulence we estimated in Sec. 2 and the 

2 If we consider the Petschek-type reconnection whose characteristic 
timescale is given by r rec ~ (L/c^) log(Rm), the magnetic field dissipation 
through reconnection might be more effective when Av/c^ < l/log(Rm) 

3 Note that the total dissipated energy increases with decreasing /3. This 
is because very small-scale corrugations develop around the current sheets in 
the case of very low plasma f3. The growth of t hese corrugations seems anal- 
ogous to the instability found by Lyubarsky (Lvubarskv 2010), and these 
corrugations slightly enhance the dissipation in the early phase. In the subse- 
quent paper of Inoue (2012), it will be shown that the interaction between 
shocks and current sheets induces instability analogous to that found by 
Lyubarsky (2010). 



study of RMI driven turbulence by Inoue et al. (2011). in 
which the velocity dispersion of turbulence is shown to be 
limited to below the sound speed. 

4. RESULTS OF 2D RRMHD SIMULATION 

Fig. Q~]shows snapshots of the structure of the perpendicular 
magnetic field, B z , at t = 0,2.5fo,7.5fo and 20fo for Av = 0.1c s , 
rj = 10~ 4 c 2 fo, k c = 3, where fo = L x /c. The turbulence stretches 
the current sheets which enhances the dissipation of magnetic 
field. 




FIG. 1 . — Snapshots of the structure of perpendicular magnetic field B- at 
t = 0,2.5t J.5t ,20t for Av = 0.1c s , r\ = 10~ 4 c 2 « , k c = 3, where t = L x /c. 

The top panel of Fig|2]shows the temporal evolution of dis- 
sipated magnetic field energy with respect to various velocity 
dispersions. The vertical axis is the dissipated energy com- 
pared to the initial magnetic field energy. In this calculation, 
we set the cut-off wave number to be k c = 3. 

The dissipated magnetic energy is calculated as AEb(& V ) = 
f B^/SirclS - f(B z ) 2 /SirdS, where Bq is the initial magnetic 
field. The top panel of Fig. 2 shows that the dissipation of 
magnetic field becomes faster as the velocity dispersion is in- 
creased. Note that the dissipation rate, however, becomes in- 
dependent of the velocity dispersion after the turbulence has 
developed and most of the magnetic energy is eventually dis- 
sipated. 

The middle panel of Fig. [2] is the dissipated magnetic en- 
ergy with respect to various resistivity, 77, with initial veloc- 
ity dispersion Av = 0.2c s and k c = 5. The solid lines are the 
numerical results with turbulence, and the dashed lines are 
those without turbulence {AEg oc 7] V 2 * 1 / 2 ). In the cases with- 
out turbulence, the dissipated magnetic energy depends on the 
square root of the resistivity, as expected. On the other hand, 
in the cases with turbulence, the dissipated magnetic energy 
seems to be nearly independent of the resistivity, in particu- 
lar, after an eddy-turnover time ~ t/to •> 9. This resistivity 
independent evolution of the dissipated magnetic energy can 
be understood as follows. 

1 . Current sheets in turbulence are gradually stretched and 
coiled by eddies, as shown schematically in Fig. [3] 
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FIG. 2. — Temporal evolution of the total dissipated magnetic field en- 
ergy. Top: the numerical calculations have been performed with respect 
to velocity dispersions Av/c s = 0.05, 0.08, 0.1, 0.3, and 0.5. Middle: 
the numerical calculations have been performed with respect to resistivity 
r]/c 2 t = 1,2,4,8, 16 X 10" 4 respectively, with Av = 0.2cj and k c = 5. Solid 
lines are the results with turbulence. Note that Av/c, = 0.5 case behaves 
somewhat differently from others in top figure. This is due to the effect of 
compression, which amplifies magnetic field in the early stage of the evolu- 
tion. Bottom: the numerical calculations have been performed with respect 
to = p gas /(B 2 /8ir) = 0.5, 1,4 respectively, with Av = 0.08c s and k c = 3. 

2. The stretching of the current sheets significantly in- 
creases their surface area, and the magnetic energy in 
an eddy dissipates quickly. 

3. This enables us to consider the effective width of the 
current sheet to be the eddy size, l tur b, indicating that 
the dissipation of magnetic energy depends only on l tur b 
instead of r\. 



5. ANALYTICAL CONSIDERATION 

In the following, we provide an analytical model of mag- 
netic field dissipation in turbulence based on the above pic- 
ture, and compare the model with our numerical simulations. 
We find that our model is in good agreement with our simula- 
tions. 
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FIG. 3. — Schematic view of the physical mechanism of the enhancement 
of magnetic energy dissipation. The bold lines show the current sheet. 
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FIG. 4. — Temporal snapshots of one-dimensional velocity power spectra 
(fP v (k)dk = J v 2 rf 3 x) in the case of Av = 0.2c,,ri/c 2 t = 1 X lO" 4 ,^ = 5. 
The spectrum quickly evolves into a Kolmogorov spectrum, P v (k) oc fc -5 ' 3 , 
well within one eddy-turnover time of the initial turbulent flow (f e ddy = 
L/Av). 

The effective width of the current sheet, l tU rb, is equivalent 
to an eddy of size corresponding to a turnover time equal to 
the elapsed time t since the induction of turbulence. In turbu- 
lence with a Kolmogorov spectrum, as realized in our simula- 
tions (see Fig. [4j, it is well known that the velocity dispersion 
of eddies of scale / is 
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where Lq is the initial energy input scale. Thus, the effective 
width of the current sheet can be written as 



turb ' 



Av 3 /V/ 2 l; 



1/2 



(7) 



The magnetic field energy in the effective current sheet is dis- 
sipated because current sheets are coiled by eddies; the dissi- 
pated magnetic energy can be evaluated as 



AE Av 



■ 1/2 



(8) 



- BlLl turh - B 2 LAv 3/2 At 3/2 /Ll 

where L is the length of the current sheet. Eq. (|8) indicates 
that the dissipated magnetic energy, AEav, is proportional to 
Av 3 ' 2 and evolves as f 3 / 2 . Fig.[5]shows the numerical results 
which confirm the above dependence. The top panel is a plot 
of the dissipated magnetic field energy AE^v with respect to 
time normalized by the eddy-turnover time. We show the re- 
sults for various Av with fixed rj and Lq. The bottom panel 
shows the dependence of the dissipation energy at t = 10fo on 
the initial velocity dispersion. Fig. [5] clearly supports our the- 
ory. 

The top panel of Fig.|5]indicates that the dissipation of mag- 
netic field is completed in approximately a few eddy-turnover 
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times. The evolutionary tracks in the top panel of Fig. [5] 
show that as the initial velocity dispersion becomes larger, 
the evolution normalized by the eddy-turnover time gradu- 
ally becomes slower. This is due to the effect of compression 
which amplifies magnetic field in the early stages of evolu- 
tion. This can be also applied to the top figure of Fig. 2 where 
Av/c.5 = 0.5 case behaves somewhat differently in the early 
phase. 



HI 

< 



0.01 




0.01 



FIG. 5. — Top: plot of the dissipation energy, AE/± V , with respect to time 
normalized by the eddy-turnover time; t e jj y = L/Av. Bottom: dependence 
of the dissipation energy at t = 10fo on the initial velocity dispersion. 



successive current sheets in a striped wind is much shorter 
than the PWN size, the dissipation can be completed suffi- 
ciently close to the inner edge of the Crab PWN. In addition, 
Mizuno et al. (2011) have recently shown that the current 
driven instability in a PWN can reduce the c-parameter from 
0.3 to 0.01. Our dissipation mechanism can also produce a 
necessary setup for their mechanism of a ~ 0.3. We also de- 
veloped a simple analytical model of magnetic field dissipa- 
tion, expressed in Eq. ||S), which agrees well with the results 
of our simulations. 

Since we assumed decaying turbulence in the postshock 
flow of a striped wind, our model can dissipate the magnetic 
field in the region where it is caught up in turbulent eddies 
with the current sheet. Hence, in order to dissipate magnetic 
field in all regions of the striped wind, turbulence should be 
driven at least on the scale of the interval between current 
sheets in a striped wind. If we consider the RMI as an ori- 
gin of turbulence, eddies generated from inhomogeneity can 
be larger than the original inhomogeneity scale. The reason 
for this is as follows: if we consider two neighboring dense 
clumps whose distance is longer than the scale of the clumps, 
the clumps can deform a shock front over their separation 
scale, which can result in eddies larger than the dense clumps 
in the downstream flow (see, e.g., Inoue et al, 2012). This may 
strengthen the applicability of our mechanism. Furthermore, 
if we account for three-dimensional e ffects, we can expect 
the induction of turbulent reconnection (lLazarian & Vis hniac 
Il999t iKowal et aL"1|2009l) . In that case, reconnection flows 
can dissipate magnetic field even if the initial driving scale 
of turbulence by the RMI is much smaller than the interval 
between the current sheets in the striped wind. Since the post- 
shock region of high-er flows still has high a Lorenz factor 
~ \/a, a single shock does not sufficiently decelerate the flow, 
and thus another shock is necessary to slow down the wind 
to match the observations of the Crab PWN. This may not be 
difficult once the turbulence converts the magnetic energy into 
thermal energy in a few eddy-turnover time. We will analyze 
this effect by performing 3D numerical simulation in our next 
papers. 



6. SUMMARY AND DISCUSSION 

In this letter, we have reported the acceleration of the mag- 
netic field dissipation by decaying turbulence, i.e. turbulence 
driven only initially by an external force in the postshock re- 
gion of a Poynting flux-dominated striped wind. Using two- 
dimensional resistive relativistic MHD simulations, we found 
that owing to the stretching of current sheets, the dissipation 
rate of magnetic field energy is independent of resistivity and 
that most of the initial magnetic energy is dissipated within 
a few eddy-turnover times, irrespective of the initial velocity 
dispersion of the turbulence. Since the interval between two 
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